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Electric Quadrupole Coupling of the Nuclear Spin 
With the Rotation of a Polar Diatomic Molecule 
in an External Electric Field 1 



By U. Fano 



Formulas are given that serve to determine the hyperfine spectrum of a rotating molecule 
in various cases, depending on the external field strength and on the rotational quantum 
number. A complete calculation is carried out for a case in which only one of the nuclei 
has a quadrupole moment, the molecule is in a rotational state (1:±1 or 0), and the inter- 
action between the field and the molecular dipole is comparable to that between the dipole 
and the nuclear spin. 



I. Introduction 

The rotation of a polar diatomic molecule in an 
electric field has been discussed theoretically by 
various authors, particularly by Brouwer 11]. 2 
Experimentally, this problem is being investigated 
at Columbia University using a modification of 
the standard molecular beam radio-frequency reso- 
nance method [2], in which the beam is subjected to 
electric instead of magnetic fields [3]. The method 
permits the energy differences between rotational 
state? of the molecule to be measured. Various 
spin-dependent interactions between the nuclei 
and the rest of the molecule may give rise to a 
hyperfine structure of the rotational spectrum. 
The electric quadrupole interaction between the 
nuclei and the molecular electrons, if present at 
all, is likely to be the major factor in this connec- 
tion and is the subject of the present theoretical 
investigation. 

Molecular rotation. — The qualitative aspects of 
molecular rotation in an electric field will be 
reviewed here, disregarding the effects that give 
rise to a hyperfine spectrum. The electronic 
system of the molecule is assumed to be in a ] 2 
state, so that it has no net angular momentum. 
The rotation of the molecule is then similar to 
that of a dumbbell whose two spheres (the two 

i Part of this work was performed at the Pupin Physics Laboratory, 
Columbia University, during the tenure of a Research Associateship of the 
Carnegie Institution of Washington in 1946. 

2 Figures in brackets indicate the literature references at the end of this 
paper. 
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atoms) carry charges of opposite sign. Several 
cases may be considered: 

(a) The electric field is so strong that the poten- 
tial energy of the molecule, when oriented along 
the field, is much greater than its rotational kinetic 
energy. Then the molecule remains oriented 
along the field, except for minor oscillations about 
this position. This case is unrealistic, however, 
as no sufficiently strong field is currently available. 

(b) The field is vanishingly weak, so that the 
molecule rotates freely about its center of mass 
with a constant, quantized, angular momentum. 
The component of this angular momentum along 
the field is then also constant and quantized. 
The positively charged part of the molecule points 
in the direction of the field about as frequently as 
in the opposite direction. Thus the molecular 
electric moment averages out, and the presence 
of the field has no effect on the rotational energy 
levels. 

(c) The field is moderately weak. The rotation 
is now appreciably perturbed by the field, and the 
molecular electric moment no longer averages out; 
an electric polarization is thus induced and each 
rotational energy level is accordingly lowered by 
an amount equal to a polarization coefficient times 
the square of the electric field. When the rotation 
has primarily the character of a precession (clock- 
wise or counterclockwise) around the direction of 
the field, the effect of the field is to force the 
positively charged part of the molecule to lean in 
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the direction of the field, and so induces a positive 
polarization. An opposite effect occurs, however, 
when the rotation brings the positively charged 
part alternately close to the direction of the field 
(downhill), and then away from it in the opposite 
direction (uphill); the positive charge moves 
rapidly while it is downhill and spends little time 
there, but it lags in the uphill position, and on the 
average it will spend more time uphill, giving rise 
to a negative polarization. 

(d) The field is intermediate or moderately 
strong. The rotation is perturbed by the field to 
a greater and greater extent. As the field in- 
creases, the positively charged part of the molecule 
is forced to lean in the direction of the field, 
whatever the type of rotation may be, and the 
polarization is always positive. 

The quantitative treatment of cases (a), (b), 
and (c) is simple, but in case (d), which obtains 
most frequently, the energy levels cannot be 
given in general as analytic functions of the electric 
field. In work as yet unpublished, W. E. Lamb, 
Jr., has developed a convenient method for the 
numerical computation of the terms for any value 
of the field strength. The rotational wave 
equation is 

[-^ 2 /2^)A-/#cos # m -W\Wm 9 ^) = 0, (1) 

where the electric field, E, is taken as a polar axis, 
(# R , <p R ) are polar coordinates of a vector/? joining 
the nuclei, I is the dipole moment of the molecule 
A its moment of inertia, A the two-dimensional 
Laplace operator (in the coordinates & R , <p R ), and 
W the energy eigenvalue. The first term in 
brackets represents the kinetic energy of rotation; 
the second represents the potential energy of the 
dipole in the electric field. A dimension! ess form 
is obtained by dividing the equation by (— % 2 /2A): 



(A+a cos &R-\-w)\p=^Q, 



(la) 



where a— IE/ (ft /2 A) is a numerical index of the 
field strength and w=W/(fr 2 /2A). The eigenfunc- 
tions can be labeled by means of two quantum 
numbers n and m; m is the magnetic quantum 
number, which indicates the rate and the direction 
of the precession of the molecule around the 
electric field; n—\m\ is the number of nodal lines, 
which characterizes the component of the motion 
along the meridians of the polar coordinates. The 
eigenvalues are functions of a. For all values of 



a, we have w n , m =w n „ m , as the energy does not 
depend on the direction of the precession. For 
a=0, the results are familiar, namely, w n , m = 
n(n+l), \f/ n> m is a spherical harmonic, and n is an 
azimuthal quantum number that measures the 
total angular momentum. For a^Cl, case (c), a 
perturbation calculation gives w nr , b ^n{n-\-\)-\- 
a 2 [(n + m) (n—m)/2n(2n—l) (2n+l) — (n+l+m) 
(n+l — m)/2(n+l) (2n+l) (2»+3)]; as expected, 
the coefficient of a 2 is positive for [m/n\ ~1 (fast 
precession), negative for |ra/w|<d. 

Quadrupole interaction. — The electrical charge 
distribution of a nucleus having a quadrupole 
moment resembles an ellipsoid of rotation rather 
than a sphere (the direction of the nuclear spin 
being the axis of symmetry of the ellipsoid). The 
value and sign of the quadrupole moment indicate 
the elongation of the ellipsoid. As the charge 
distribution of a diatomic molecule is not spheri- 
cally symmetrical about any of its nuclei, the 
electrostatic energy of the nucleus within the 
molecule will depend upon the orientation of the 
nuclear spin with respect to the axis of symmetry 
of the molecule (i. e., the line joining the nuclei). 
If the spin be directed along the axis and then 
swung into the opposite direction and back into 
the initial one, the interaction energy will go 
through two maxima and two minima. 

In the absence of this quadrupole interaction 
and of magnetic interaction, the nuclear spins 
would maintain a constant orientation in space 
while the polar molecule rotates in an electric 
field, and the orientation of each nuclear spin 
with respect to the axis of the rotating molecule 
would undergo periodic variations. In the pres- 
ence of a quadrupole moment there will then be 
periodic variations of the interaction, and hence 
torques that affect both the rotation of the mole- 
cule and the orientation of the nuclear spins. This 
effect will be particularly conspicuous in case (b) 
described above, as the orientation of the angular 
momentum of the molecular rotation can then be 
shifted by a torque without perfoiming any work. 
This angular momentum will then precess freely 
together with the nuclear spins, and only the total 
angular momentum, i. e., the vector sum of the 
molecular momentum and of the nuclear spins, 
will remain constant. When, on the other hand, 
considerable work must be performed to upset 
the molecular rotation, the quadrupole interaction 
will have no significant effect on the motion of the 
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molecular constituents; the energy of the system 
will then simply include a contribution from the 
quadrupole interactions averaged over the varying 
orientation of the spins with respect to the mole- 
cular axis. This contribution will depend on the 
state of rotational motion of the molecule and on 
the (constant) orientation of the nuclear spins 
with respect to the electric field. This situation 
prevails generally in cases (a) and (d) described 
above, i. e., for moderately to very strong fields; 
from this standpoint, case (c) represents an 
intermediate situation. 

It should still be considered that no work is 
ever required for the simple operation of reversing 
the direction of the precession of the molecule 
around the electric field, provided the speed of 
precession is unchanged. Owing to the quadrupole 
interaction, and for any electric field b1 length, the 
uniform precession might therefore conceivably 
be turned into a precession whose direction is 
periodically reversed; the conservation of angular 
momentum would be preserved by corresponding 
variations of the precession of the nuclear spins. 
Quantum mechanically, there would thus be 
"standing waves of precession/ 7 the number of 
nodes of the standing wave depending on the 
speed of precession. However, it will be shown 
further below that only the standing wave with 
two nodes actually arises in this way. This is a 
direct consequence of the particular dependence 
of the quadrupole interaction upon the mutual 
orientation of the molecule and of the nuclear 
spin; as already stated this interaction goes 
through two maxima and minima when a cyclic 
change of orientation is performed. 

The quantitative treatment of our problem 
will accordingly proceed as follows. As shown by 
Casimir, [4], the interaction energy between a 
nuclear quadrupole moment and the electronic 
charge of an atom or molecule is of the order of 
e 2 qQ/4: (q and Q are the usual symbols for a 
molecular and a nuclear constant whose product 
is <1 cm~ l ). Hence, taking into account the 
possible values of the pertinent constants, it 
seems safe to assume for the dimensionless 
parameter characterizing this interaction: 



\0\=\e\QIHPl2A)\<lO- 



(2) 



This value must be related to the fact that, in gen- 
eral, the separation of different unperturbed rota- 



tional levels of the molecule is 3 |t*V W '— w nw |>l. 
Therefore, in general, the quadrupole interaction 
can be treated as a small perturbation of each 
rotational term, and each term will preserve its 
quantum numbers n, m. This procedure will not 
be adequate in the cases discussed above, in which 
the quadrupole interaction has a great effect on 
the rotation of the molecule, namely, for states 
(n, m) which are directly coupled by the quad- 
rupole interaction to other states (n', m') such 
that |w»'m'— w nm |<|0|. Direct quadrupole cou- 
pling means that the quadrupole interaction energy 
H Q has a matrix element (n'm'\H Q \nm)^0; this 
occurs only when \m— ra'|:<2, as it will be shown 
in the next section. These two conditions can be 
fulfilled simultaneously in two cases only: 

(1) n=n'>0, m=—m'=±\ (all values of a), 
in which the molecule processes back and forth, 
i. e., clock- and counterclockwise; 

(2) n=n'>0, c* 2 <|/3| (|ro-ro'|^2), in which the 
field is weak and the rotation is strongly perturbed. 

The latter case becomes simpler when a 2 <c|/3| or 
even a 2 =0. In fact, the zero-field case involves 
a simple vector coupling of the rotational angular 
momentum J (J=n) with the nuclear spins I x 
and / 2 . In the case of near-zero-field, this vector 
coupling is preserved and the effect of the electric 
field can be treated as a small perturbation, some- 
what analogous to the anomalous Zeeman effect 
produced by a weak magnetic field. 

The general expression for the quadrupole 
interaction matrix is given in the next section; the 
application to the different cases will then follow. 

II. Matrix Elements of Quadrupole Inter- 
action Between Molecular Rotation and 
a Spinning Nucleus 

(a) The electrostatic interaction — e 2 /\r— p\, 
between a nuclear proton at a position p with 
respect to the center of the nucleus and a molecular 
electron at a position r with respect to the same 
center, can be expanded into a series of Legendre 
spherical functions of the angle y rp between the 
vectors r and p. The first term of this series 

3 The magnetic spin-spin interaction between the nuclei and the inter- 
action between the nuclear magnetic moments and the rotation-induced 
magnetic moment of the molecule are generally weaker than the electric 
quadrupole interaction; in our scale they may be characterized by a parameter 
»?<10- 5 . These interactions will be disregarded in the following therefore 
we must assume 
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represents the Coulomb attraction of the electron 
toward the center of the nucleus, the second term 
has no average effect, and the third, namely: 

H Q =-e 2 (p 2 /r*)P 2 (C0S 7pr ) = 

-e>7r 3 )(|cos 2 Tpr -^ (3) 

is called the quadrupole component of the inter- 
action. This component must be integrated over 
the density of protons throughout the nucleus 
and the density of electrons throughout the 
molecule. The integration over the electron 
density is carried out in polar coordinates having 
as axis a vector R from the center of the molecule 
to the center of the nucleus. Using the addition 
theorem of spherical harmonics and taking into 
account the fact that the 2-type electronic state 
has axial symmetry, we find: 

<(l/r 3 )P 2 (cOS 7 pr)>r = 
((l/r 3 )P 2 (cOS 7*r))rP 2 (cOS 7/j p ), 

where the symbol () r indicates integration 
with respect to the electron density. The first 
factor on the right-hand side is a molecular char- 
acteristic which is indicated [4] by q/2. Using 
again the addition theorem of spherical harmonics, 

P 2 (cosy Rp ) = 

(47r/5)S M (-l)MF 2 ,_ M (^ p , <p P )Y 2t t(# M , <Pm), 
we have, 

(H Q )r = 

(^V2)2 m (-1MpW5) 1/2 F 2 , _ M (tf p , <p p )] 

X[(4tt/5) 1 / 2 F 2 , m (^, <p R )\. (3a) 

(The normalization of the spherical harmonics is 
suchthat//|r ? , M | 2 sin^#^-l,r^=(-l)^r^_ /i . 
(b) All the matrix elements of the first factor in 
square brackets, i. e., the integrals of this factor 
over the proton density, associated with pairs of 
nuclear states with the same spin I and magnetic 
quantum numbers m 7 , m/, can be given in terms 
of just one such matrix element. Thus: 

(m/b 2 (4T/5) 1 / 2 F 2 , _ M (# P , < Pp )\m I ) = 
(7b 2 (47r/5) 1/2 F 2 , (0 P , *)\I)(P#IDZ)*mS. -thh 
The first factor on the right-hand side is a nuclear 
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characteristic which is indicated [4] as Q/2 (i. e., 
half of the quadrupole moment). In the formula 
8 it k =l or 0, depending on whether i=k or i?±k. 
The coefficients D)'™' may be found by means of 
group theory, 4 and are given in table 1. As 
Z>JJ = — 7(27— 1 )/ (27— 1 ) (27+3) we take: 

d(j)t= ~(2j-l)(2i+3)J%'; 

the coefficients d are given in table 2. Thus: 

(m 7 -Mb 2 (47r/5) 1/2 r 2 , »(#p, <P P )\™r) = 

Qd(I)% r -»/2l(2l-l). (4) 

Table 1. — Coefficients D];™ for: j'—j = y ±2; 





m' — m = 0, ±1 


, ±2" 




D) 


:=-[3m2-j(j+l)]/(2j-l)(2/+3) 






V) 


:: +1 = (2m+l)[(3/2)(j+m-fl)(j-r«)l l 


4 




D) 


^ +2 =-[(3/2)0'+w+2)a+7»+l)(j-i 


n)(j—m- 


-U]i/2/(2;-l)(2;+3) 


X>^- m = (3/2)[(i+7»+2)(i+m+l)(j-m- 


f2)0"~» 


+ l)/(2j+l)(2/+3)2 




(2j+5)]'/2 






D \ 


h i '" +1 = [(3/2) (j+m+3) (jj-m+2) (j+1 

(2j+5)]'/2 


«+l)0- 


w+l)/(2j+l)(2j+3)2 


»\ 


^• w,42 = [(3/8)(i+w+4)(j+w+3)(;+m+2)(j+ 


w+l)/(2j+l)(2j+3)2 




(2/+5)] 1 ' 2 







° All other coefficients are obtained by means of the formulas 

Table 2.— Coefficients d(j)Z' = - {2j- 1) (2j+B)D% 



<c 2 


= KW)(J+m+2)U+m+l)(j- 


-m)(j—m- 


-l)) in 




<c +1 


= -(2m+\)\&2)(j+m+V(J- 


-m)\ v2 






*z 


= 3m2-J(j+l) 








^r 1 


= (2w-i)[(3/2)o-m-fi)a+w)] l/2 =rf:;: 


H 




d m-2 


= H3/2)(j-m+2)(j-m+l)(j-\-m)(j+m- 


-ni ,/2 = 


=d:T 2 



(c) The matrix elements of the factor in the 
second bracket of eq 3a, pertaining to the rotational 
states of the molecule (n, m) and (n' , m') is the 
average of this factor over the probability distri- 
bution ^*/ t m ^ n , m and can be indicated as 

* These coefficients are closely related to those indicated by E. Wigner, 
Gruppentheorie (Vieweg, Braunschweig 1931), p. 206, as £$?„. For j and f 
integer D]'™' = (4x/5)WffY*j'- m >(&<p) Y 2 . «'-«(*, ?) Y,- , m (»<p)simddd<p. The 
formula given means essentially that the ratio of matrix elements for 
different values of (m/, mi') is the same as though the angular distribu- 
tion of the proton density in each nuclear state (7, mi) were described by 
the spherical harmonic Fj, m/ (suitably generalized for half-integer 7). 

Journal of Research 



n', m'|(4,r/5) 1/2 F 2 „(0«, <p K )\n, r») = (4,r/5)"W>;. 

Expanding the rotational wave functions ^ n , m 
into spherical harmonics: 

* n , m (6, <p)=2 3 C%"Yj, m (6, «,), 

where the coefficients Cj m depend on the parameter 
a, we have: 

T>n',m' y (TlJ + 2,m'f > in'm'\ TUmTln'm' I 

■W:l' m 'Cl^)CT. (5a) 

For n= n' ', m=m' (which is the only case 
involved in a first-order perturbation), we have: 

Rl:2=R(n y m) = (3/2)(cos 2 d) R -l/2. 
There is also a simpler form for 

R n n :V = S(n) = (3/8y< 2 (sm*d R ). 

The coefficients R must be obtained by numeri- 
cal computation except when a is very large or 
very small. For a<^Cl: 

Cj' m =8n t j+ot&n*hjf(n, ?i±l) + terms in a 2 , etc., [I] 

where / is a function of n and n±l; hence eq 5a 
gives 

Rlz'=Dl-% + terms in a 2 , etc 

For cCM, l^n./wl^O, except for # R <Cl (i. e., the 
molecule points approximately in the direction of 
the field). Hence, according to eq 5, lim #£'£*'= 

f>n,n f 8m,m'- However, the convergence to this limit 
is slow. For example, R(n, m) ean be expanded 
in powers of a~ l/2 as follows: 

R(n, m) = 

l~3(27i-m+l)(2a)- 1/2 + terms in a" 1 , etc. 

(d) Combining eq 3a, 4, and 5, the matrix 
element of 7/ C; , pertaining to states with the same 
total magnetic quantum number M=m I J rm = 
m f I -\-m / , is found to be 

(n/, m'\H Q \n, m) = 

- (« 2 #) (- i) m/ - m R^d(iy A f r-Z'/i(2i-- 1) - 

- (#72^1) (- l)^"'jr%T7J(7)iJi:JV (6) 



,/ (0h<p;t)Y 2jl (d Ii , <p:t)f„, m (OR, <pr) sin d K dd R d<p . 



where 



(5) 



^ = |8//(2/-l)=[^(2/4/(2/~l)]/(fe 2 /2^). 
III. The General Case: a 2 »|j8|, |ra|^l 

Here it is sufficient to treat the quadrupole 
coupling as a first order perturbation of each 
rotational term. 

Only one of the nuclei has a quadrupole moment. — 
Each term of the complete system is characterized 
by the three quantum numbers n, m, and 
m 7 =M— m. The term energy is then given, to 
the first order in £, by the sum of the rotational 
term and of the matrix element (6) with n=n' 
and m=m': 



W n , m , m = (h 2 /2A) [w n , m - p'd(I)ZiR n n%} = 
(P/2A){w n , m -(3'[Sin 2 I -I(I+l)mn,m)} = 



W n 



r YV n , 



t "' n ,— r 



(7) 



Each rotational energy level of energy w n ^ m \ splits 
thus, in this approximation, into /+§ levels if I 
is a half-integer, into 7+1 levels if 7 is an integer. 
The absolute values of the intervals between the 
levels of such a multiplet are functions of the 
electric field (through 7?), but their ratios are 
independent of it . 

Both nuclei have a quadrupole moment.— Each 
term is now characterized by four quantum num- 
bers n, m, m IU m I2 . The term energy is, again to 
the first order in /?: 



W\ 



n , ± m, =*= m ji , ± mil 



(h 2 l2A){w n , m -^[Sm 2 I1 -I l (I 1 +l)}R(n ) m)- 
e 2 [3m 2 12 -I 2 (I 2 +l)]R(n,m)} (8) 

IV. The Special Case: a 2 »|/3|,|ra| = l 

Only one oj the nuclei has a quadrupole moment. — 
The energy of terms with m=l, ra 7 — (7— 1) or 7, 
M=I or (7+1), or with m.= — 1, m r = — (7— 1) 
or —7, M=—I or —(7+1), is still given by eq 7. 
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For other terms, with \M\<^I, m is no longer a 
good quantum number, but \m\ is. The term 
energy e(fi 2 /2A) is then obtained, to the first order 
in 0, by solving the secular equation, 



=0, 



w nl +(n,l\H Q \n,l) — e (n,l\H Q \n } —l) 

(n—l\H Q \n,\) Wn^t+in,— l)\H Q \n,— 1) — e 
with R(n,l)=R(n, — l), w n .i=w»,_i. Its solution is 



[(3/2) (I+M) (I+M+ 1) (I-M) (I-M+ l)S(n) 2 /R(n f l) 2 + 36M 2 ] 1 ' 2 } . 



(9) 



The original w nl rotational level is thus split into 
27+ 1 different levels, if 7 a half integer, and into 
2(7+1), if / is an integer. 

Both nuclei have a quadrupole moment. — This 
case is similar to the preceding one, but more com- 
plicated. There are three levels with total mag- 
netic quantum number \M\=I 1 +I 2 or 7i+7 2 +l, 
for which the energy can be calculated by means 
of eq 8. These energies are indicated by the 
following symbols: 

vv n .1 ,I 1 ,I 2 ~ ^ n ~ l - J l ~ J 2 » 

W n ,i ,/j j 2 ~i : = ry n> _i < _/ lt _7 2 -|_i ; 

yyn,i,i 1 -i,i 2 ==: W ni -i-i x +i-j 2 , 

For each group of terms with |MK7i + 7 2 one 
must solve two secular equations that in general 
will have more than two rows and columns. The 
rows and columns correspond to zero-order states 
characterized by particular values of (m, m n , m I2 ) 
with |m| = l, m+m n -\-m I2 =M. The zero-order 
states belonging to each secular equation can be 
arranged in a succession, which (for AT>0, 7 2 >7i) 
can be of either type: 

(1,7,, M-1,-1), (-l,I h M-1, + 1), 

(l,/ 1 -2, M-7^-1), (-1, 1,-2, M-1,+3) . . . 



or 

(l,/i-l, M-I,), (-1,/i-l, M-It+2), 
(1,^—3, M-I x +2), (-1, 1,-3, M-1,+4) . . . 

The diagonal elements of the secular equation are 
of the type: 



/mi2i 



w nl -R{n,l)W\d{I l )Z\\+^ 2 d(I 2 )l 
the elements next to the diagonal are of the type: 

-SinW.dd,)^ 2 or -S(n)t3' 2 d(I 2 )Zr 2 
as the case may be; all other elements vanish. 

V. The Special Case: a<l, n>0 5 

Only one nucleus has a quadrupole moment. — 
For a<Cl, or, more accurately, for: 

a 2 /2n(n+l)^[al(w n + ltm -w n , m )] 2 <\p\<^l 

the difference R^—D^^a 2 can be disregarded, 
since R is multiplied by f$ in the expression of the 
matrix elements of quadrupole interaction. 

Replacing then R with D and using the approxi- 
mate expression for w n>m valid for small a 6 , the 
matrix element of the complete energy operator 
pertaining to states with the same values of n, 7, 
and M==m+m r is 



{m'\H\m)^<p/2A){n(n+l)t mtM .-- 

^(-l)^- m l):I(l(^:;7/(2/-l)-5 w ,a 2 /2/i(n+l)^ 
terms in a 4 , a 2 fi, or 2 , etc. 



The term energy e(ii?/2A) is then obtained, to the 
first order in 0, by solving the secular equation: 



\(m'\H\m) — e(ii 2 /2A)5 r , 



= 0. 



The number of rows and columns in the determi- 
nant, and hence the degree of the equation, is 

8 For a<^l, w=0 we have BH^a 2 ; hence the quadrupole interaction can 
be disregarded. 
6 This expression can be given in terms of the coefficient D™* provided w>0. 



^(2n+l) or ^(27+1), whichever is the smaller, 
since m and ra' take all the values that do not 
exceed either (M+I) or n and are not less than 
either (M—l) or — n. The term energy can also 
be indicated as: 



(ft 2 /2A)[(n(n+l) + p"y)], 



(10) 



where P" = (3/(2n-l)(2n+ 3)7(27-1) and y is a 
solution of the secular equation 



220 



Journal o£ Research 



i (- l) m - m/ d(?Om (d(I)%-2'+d m > , m x) -d m > , m y\=0 

(10a) 

Here x = -[7(2/-l)/2n(n+l)](c* 2 /0). The coeffi- 
cients d are given in table 2, hence the further 
steps in calculating y for any given value of x 
(i. e., of the field strength) are* straightforward 
even though they may be laborious. The particu- 
lar case for 7i= 1 is worked out in the appendix, 



which also contains some comments about the 
methods of solving the secular equation. 

Both nuclei have quadrupole moments. — The cal- 
culation is similar to that of the preceding case 
but is more complicated. The energy of each 
term can be expressed as: 

(P/2A)[n(n+l)+y-$l(2n-l)(2n+3)]. (11) 
Where y is a solution of the secular equation: 






= 0. 



(11a) 



Here x= — a 2 /p2n(n J rl) and 0, b x and b 2 are so 
chosen that: 

ff5i=ft//i(2/ 1 -l) l 5fc=A//,(2J 1 -l) l & 1 ~6 2 ~l. 

Finally, |ro|,|ro'|^n; \mnUm'n\^Iu\rn I2 \,\m' I2 \^l2; 

7n-{-m n -{-m I 2 = m / + m'ji + m'j2=M. 

VI. The Special Case: a 2 <|/3| 

Zero field. Only one nucleus has a quadrupole 
moment. — In this case the total angular momen- 
tum F=J+I (J=n), i. e., the sum of the rota- 
tional angular momentum and of the spin of the 
nucleus having a quadrupole moment, is a con- 
stant of the motion of the system. J and / are 
not constant but precess freely around F on 
account of the quadrupole interaction. All the 
states having the same value of F and different 
values of the total magnetic quantum number 
M=m J rm r have the same energy. As shown by 
Casimir [4] 7 the energy of each term is propor- 
tional to an eigenvalue of the quadrupole coupling 
operator: 



| (/•!)( J-D+|(/.i)- 



1 



(/•/)(/•/). 



The proportionality factor can be determined by 
comparison with the solution obtained in section V 
for the case M=F=n J rI, a 2 =0, which is trivial. 

The energy levels for n-\- I^F>:\n — 1\ are 

(P/2A {n(n+l)+fflC(C+l)-4n(n+l)I(I+l)] 
/2(2n-l)(2w + 3)J(2l — 1)}, (12) 

with 

C=2(J.I)=F(F+l)-n(n+l)-I(I+l). 

Near-zero field. Only one nucleus has a quad- 
rupole moment.- — The weak electric field has the 
effect of splitting the levels with equal F and 
different M. The energy is calculated by con- 
sidering the term proportional to a 2 in eq 10, 
that is, 

a 2 DlZ/2n(n+l) ; 

as a small perturbation. This means simply 
averaging the value of this term over the partici- 
pation of each rotational state (n, m) in the com- 
plete state (n, F, M). The following term must 
therefore be added to eq 12: 



(n 2 /2A)[a 2 J2n(n + l)]: 



(fHnl) 



^ m .M-J 2 D:^-a 2 (P/2A)^M 2 -F(F+l)]{3D(D-l)-4F(F+l)n(n+l)]/ 
2n(n+l)(2n-l)(2n+3)2F(F+l)(2F-l)(2F+S), 



(13) 



_„ is the coefficient of the eigenf unc- 
tion (n, m; I, M—m) in the linear combination 
which constitutes the eigenfunction (n, F, M) [5], 
and 

D=2(J.F)=F(F+l)+n(ri+l)-J(I+l). 



i H. B. O. Casimir, Archives du Mus'e Teyler [III] VIII, 201 (1936). The 
application to a molecule in zero-field has been made by B. T. Feld and 
W. E. Lamb, Jr. (Phys. Rev. 67, 15 (1945)). The definition of q by Feld and 
Lamb is, however, different from that introduced in this paper; this entails 
a difference in the multiplicative factor of some formulas. 



In the theory of the anomalous Zeeman effect 
[6] for a case of (LS) coupling, which is somewhat 
analagous to the present calculation, one has to 
calculate the average value of the quantum 
numbers m L and m s , namely, ^ mL (Cf t ^ lLy M _ mL ) 2 m Ly 
etc. It is shown by direct methods, however, 
that the results are M(L-J)/J(J+1) and 
M(S-J)/J{J+1). In our case the quantity D™ 
is a matrix element of quadrupole interaction, 
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while m L is a matrix element of dipole interaction. 
It might therefore be expected that eq 13 could 
be expressed directly by means of quadrupole 
matrix elements. In fact eq 13 involves the 
expression D F F % and the matrix element of the 
operator: 

| (/.F)(/.F)-f (/-F)-| (J-J)(F-F), 

but the latter differs from the usual quadrupole 
coupling operator by the sign of the second term. 
Both nuclei have quadrupole moments. — This case 
is again similar to the preceding ones but more 
complicated. Equation 13 should be extended to 
include two terms, one involving the coupling of 
/ and I u and the other that of J and / 2 . However, 
neither «/+/i nor «/+/ 2 is a constant of the 
motion, but only the sum F=J J rI 1 J rI 2 . In gen- 
eral, there will be more than one term with the 
same values of n and F. The characterization of 
each term and the evaluation of (J-Ii) and (J-I 2 ), 
and hence of the energy levels, depend upon the 
ratio ft /ft and require the application of methods 
similar to those used in the theory of complex 
atomic spectra involving a number of terms with 
the same angular momentum [6, p. 233]. The 
same holds for the extension of equation 13. 



=0, 
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VIII. Appendix 

The solution of the secular equation 10a for n = l is 
obtained as follows. The total azimuthal quantum 
number M assumes all the values |il7|:<7-f- 1. For 
|il7| = 7+l the secular equation has only one row and 
column; its trivial solution is: 



y = I(2I-l)+x, for |M| = 7+1. (14a) 

For | M \ = I, the equation is quadratic: 

[(7-3) (27- l)+x]-y 3(2/-l)/V2 
3(27- l)P/2 -2[7(27- 1) +x]-y 
and its solution is 

y=~ !{(/+3)(2l-i)+3=fc 

3([(7-l)(27-l)+x] 2 -f-4[7(27-l)p)i/2|, for \M\ = I. (14b) 

For |A/|<7, the equation is cubic, having three rows and 
columns. Its solution can be indicated as 

y = K l + 2K 2 v* cos [(a + 2p*-)/3], for |M|<7, (14c) 

where p = 0, 1, 2, cos a=-Ki/K<?#. K u K 2 , 7T 3 are 

polynomials of first, second, and third degree in the 
elements of the determinant form of the secular equation, 
which can be given in terms of 7, M and x: 

K x = 2 

K 2 = [4P(I+ l) 2 + 7(7+ i)+i] + 2[3il!2-J(/+ 1) + l]x+x* 
K 3 =[8P(I+ l) 3 -247 2 (7 + l)2-h (3/2)7(7+ 1)]- 

3{3ikf2[2/(/+l)-7/2]-J/(/+l)[4/(/+l)-r-l]-M^+ 
3[32kf 2 -I(J+l)-f lja^+s 3 . (15) 

Figure 1 shows the plot of y as a function of x for 7 = 7/2, 
corresponding to the Cs nucleus in the CsF molecule 
(F has no quadrupole moment), and for all values of M. 
The position of the energy levels for x greater than two or 
three times 7(27—1) is already given with good accuracy 
by eq 7 and 9, which were calculated on the assumption 
of a 2 ?>>\8\. 

The derivation of the eq 15 is a purely algebraic but 
rather tedious and laborious process. This is probably 
true in general of the process of solving secular equations 
in which the matrix element are still functions of various 
pertinent parameters. It is not just the solution of the 
equation that is tedious, but especially the preliminary 
process of transforming the equation from its determinant 
form to a polynomial form. The labor involved increases 
rapidly with the degree of the equation. Therefore, when- 
ever it is necessary to obtain a solution only for a few 
values of the parameters 7 and M and x, it is advisable to 
introduce these values in the matrix elements themselves, 
and to carry out the rest of the work numerically. Other- 
wise it may be advisable to consider other ways of reducing 
the equation to a polynomial form. For example, the 
coefficients of an algebraic equation can be expressed in 
terms of the sums of powers of its unknown solutions y im 
In the example above, Kt=Si/B, # 2 =(ft-£?/3)/6, 
K s =((S,-S 1 S^/6)-S 3 1 /27, where S„=2^. The quan- 
tity S n is the trace of the matrix of the n-th power of the 
energy operator; it may well be simpler and safer from 
mistakes to calculate the necessary S n than to expand the 
determinant form of the secular equation. 
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Figure 1. — P/o/. o/ /fie energy levels W n versus the square of the electric field, E 2 , for n=l, 1=7/2, \M\.< /+ 1=9/2; rwir- 

dinates are in dimensionless units: 

* = -21 (lE)ye 2 qQ(hy2A), y = 42QW/e*qQ. 
— , M=±l/2; .... Jlf=±3/2; Af=±5/2: _._._., A/=±7/?; ++++, M-=fc9/2. 

Washington, March 31, 1947. 
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